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We calculate the skewed quark distributions(SQDs) of the 
pion in the light-front quark model, and discuss the calcula- 
tion of the nonvalence contribution to the SQDs in this model. 
The frame-independence of our model calculation is guaran- 
teed by the constraint of the sum rule between the SQDs and 
form factor. Our numerical results show large nonvalence con- 
tributions to the SQDs at small momentum transfer region as 
the skewedness increases. 

PACS number(s): 13.40.Gp, 13.60.Fz, 12.39.Ki 



I. INTRODUCTION 

Recently there has been a great interest in the ofF- 
forward (or nonforward, off-diagonal) parton distribu- 
tion functions defined as non-diagonal hadronic matrix 
elements of bi-local products of the light-front quark 
and gluon field operators @,|,|,|j5|§,|,|,|jl0) . The 
off-forward parton distribution functions, the so called 
"skewed parton distributions (SPDs)", are the general- 
ization of the ordinary (forward) distribution functions. 
A well-known and practical example of SPDs as a nonper- 
turbative information entering the light-front dominated 
hard scattering processes is the deeply virtual Compton 
scattering (DVCS) j*p — *■ 7p for large initial photon vir- 
tuality and small t region, which can be factorized 
into a hard photon-parton and a skewed parton distribu- 
tion |l|,|,|. 

In the present work we formulate the pion form factor 
in terms of the SPDs. Since the usual local photon vertex 
in the pion form factor analysis is replaced by a nonlocal 
operator of the SPDs, one can explore new physics. The 
physical interpretation of the SPDs becomes clear in the 
light-front frame with the light-front gauge yl+ = 0. In 
the light-front coordinates, the SPDs are in general func- 
tions of the longitudinal momentum fraction variable x, 
the skewedness parameter ^ — {P — P')^/P^ measur- 
ing asymmetry between initial (P) and final {P') hadron 
state momenta, and the squared momentum transfer t. A 
simple physical interpretation of the SPDs, in light-front 
quantization, is that they provide a link between the or- 
dinary parton distributions of hadrons and the hadronic 
form factors, i.e. the ordinary parton distributions arc 
forward (^ = and t = 0) limits of the SPDs and the 
form factors are given by moments of them. 

Due to this dual role of SPDs, they are closely re- 
lated to form factors with the only difference between 
SPDs and form factors being that the momentum of the 



"probed quark" in SPDs is not integrated over but rather 
kept fixed at momentum fraction x. For example, in 
studying the light-front wave functions of hadrons, the 
overlap representation of the light-front wave functions 
of hadrons for the form factors in spacelike {q^ < 0) re- 
gion can be obtained if one uses J^{= J° -I- J^) and 
the Drell-Yan-West (g"*" — 0) frame where only parton- 
number- conserving valence Pock state contribution is 
needed. The successful phenomenological calculations 
of form factors in spacelike region can be found in the 
light-front quark model(LFQM) |ll|,|l|,|l|,0 . On the 
other hand, form factors in timelike (q^^> 0) region such 
as the weak form factors for exclusive semileptonic de- 
cays require g"*" > frame, which in turn require parton- 
number-changing nonvalence Fock state contribution as 
well as the valence one [ p^Hl^ . Similarly, while the or- 
dinary parton distributions (analogous to = limit 
of form factor) can be represented in terms of squared 
light-front wave function of a hadron, one cannot avoid 
nonvalence contributions to SPDs since they always in- 
volve non-zero ^ corresponding to g"*" 7^ in timelike 
form factor calculations. In recent papers 

10, the non- 
valence contribution to the SPDs has been rewritten in 
terms of light-front wave functions with different par- 
ton configurations. However, the representation given 
by |^,|l^ requires to find all the higher Fock-state wave 
functions while there has been relatively little progress in 
computing the basic wave functions of hadrons from first 
principles. Our approach provides an alternative way of 
handling the nonvalence contribution which is more suit- 
able for the constituent quark model(CQM) specific to 
the low momentum transfer processes. We use the light- 
front Bethe-Salpeter (B-S) formalism for the SPDs as an 
extension of our treatment of the pion form factor [ p^ . 
Although the present work has many features similar to 
that of the Fock state expansion, we utilize the close re- 
lation in the B-S formalism between what is interpreted 
as the valence and nonvalence wave functions in the Fock 
state approach. 

In an effort to apply light-front wave function based 
phenomenology to form factors in timelike exclusive pro- 
cesses, we have presented in [|l^ an effective treatment of 
handling the nonvalence contribution to the weak form 
factors, based on the B-S formalism, and obtained rea- 
sonably good numerical results for the processes in the 
small momentum transfer region. The main purpose of 
the present work is to apply the effective method pre- 
sented in MM of handling the nonvalenece contribution 



1 



to the SPDs of pion at small momentum transfer region 
in LFQM. The paper is organized as follows. In Section 
II, we briefly introduce the necessary kinematics in which 
we follow the notation employed by Radyushkin |3|. In 
Section III, we represent the SPDs of the pion in terms 
of light-front vertex functions, starting from the covari- 
ant Bethe-Salpeter model of (3 + l)-dimcnsional fermion 
field theory. The nonvalence part of SPDs is expressed 
in terms of light-front vertex functions of a hadron and a 
gauge boson. The link operator connecting {n — l)-body 
to (n-l- l)-body in a Fock state representation is obtained 
by an analytic continuation of the usual B-S amplitude. 
We also show that the complicated (n -I- 2)-body energy 
denominators are absent through the calculation of the 
light-front time-ordered diagrams. Of particular interest, 
the instantaneous contribution of the quark propagator 
to the nonvalence diagram is seperated from the on-shell 
propagating part. In Section IV, we replace the light- 
front vertex functions obtained from Sec. Ill with our 
LFQM wave function [|l2| and show our numerical results 
for the SPDs of the pion at small momentum transfer re- 
gion. We also show that the frame-independence of our 
model is guaranteed by the sum rule between the SPDs 
and form factor of the pion. Conclusions follow in Section 
V. 



Defining the four momentum transfer A = P — P' , one 
has 



P = 
P' = 



nf2 I A 2 



and 



A = P- P' 



a -OP' 



^ A2 + A2 



(2) 



(3) 



where M is the pion mass and ^ = A+/P+ is the 
skewedness parameter describing the asymmetry in plus 
momentum. The squared momentum transfer then reads 



t 



2PA 



(4) 



Since A^ > 0, t has a minimum value — tmin = 
^/ (1 — ^) at given ^. As shown in Fig. the parton 
emitted by the pion has the momentum fc, and the one 
absorbed has the momentum k' . 

As in the case of spacelike form factors, we choose a 
frame where the incident spacelike photon carries ^ 0: 



II. KINEMATICS 



hard 




soft 



FIG. 1. Handbag diagram contributing dominantly to 
Compton scattering in the deeply virtual region. The lower 
soft part consists of a hadronic matrix element which is 
parametrized in the form of skewed parton distribution func- 
tions. 

We begin with the kinematics of the virtual Compton 
scattering (see Fig. ||) of the pion 

7*(g) + ^(P)-7(9') + 4n, (1) 

where the initial (final) hadron state is characterized 
by the momentum P (P') and the incoming spacelike 
virtual and outgoing real photon momenta by q and 
<?', respectively. We shall use the component nota- 
tion V = {V~^ ,V~ ,Y±) and our metric is specified by 
F± = {V° ± V^) and = V+V- - V^. 



0, 



(q± + Aj 



^p+ 

(q± + A±)^ 



{1-0P+ 



(5) 



In deeply virtual Compton scattering (DVCS) where 
= —q^ is large compared to the mass M and —t, 
one obtains 



2P-q 







(6) 



i.e. ^ plays the role of the Bjorken variable in DVCS. For 
a fixed value of — t, the allowed range of ^ is given by 



< 6 < 




(7) 



III. SKEWED QUARK DISTRIBUTION OF THE 
PION 

Analogous to the pion electromagnetic (EM) form fac- 
tor calculation 

j+(o) ^ {p'\mi^m\p) = F.mp + pt, (s) 

we define the skewed quark distributions(SQDs) 
J-Tt{£,, X, t) of a pion by 
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(9) 



where z = (z"*", z_l) in a light-front representation. 
Note that the path-ordered exponential of the gauge field, 
V exp[i J z^Afj], required by gauge invariance in Eq. (|^) 
does not appear in the light-front gauge yl+ = 0. As one 
can see from Eqs. (||) and (^), the involves one less 
integration than the form factor Fj^ due to nonlocality 
of the current matrix element. The SQDs display char- 
acteristics of the ordinary(forward) quark distribution in 
the limit of ^ ^ and t — > 0, on the other hand, the first 
moment of the SQDs is related to the form factor by the 
following sum rules : 



Jo 



(10) 



where :r^(f,a;,i) = e„.f^(^,x,i) - edT^{^,x,t) and we 
assume isospin symmetry (m„ ~ mj) so that x, t) = 

X, t). Note that Eq. ( p^ is independent of ^, which 
provides important constraints on any model calculation 
of the SQDs. In general, th e poly nomiality conditions for 
the moments of the SQDs defined by 



Jo 



(11) 



require that the highest power of ^ in the polynomial ex- 
pression of Fn{^,t) should not be larger than n. These 
polynomiality conditions are fundamental properties of 
the SQDs which follow from the Lorentz invariance. In 
case of a spin- 1/2 composite system [§|J8), where the he- 
licity non-flip(_ff ) and helicity flip(£') SQDs are involved, 
the second(n = 2) moment of each SQD yields the ^ 
dependent form factors of the energy-momentum ten- 
sor, although the sum of the moments produces the ^- 
independent form factors of the energy-momentum ten- 
sor. When the sum rule(n = 2) is extrapolated to 
— t = it provides the information on the total 
quark(i.e. quark orbital) contribution to the nucleon 
spin. For the spin-0 composite system like the pion, 
the situation is quite different because only one SQD 
X, t) exists. We discuss our numerical results of 
Pn{^, t){n — 1, 2, 3) for the pion in the next section (Sec- 
tion IV). 



FIG. 2. Diagrams for SQDs in different kinematic regions 
for the case ^ > 0: The covariant diagram (a) corresponds to 
the sum of the LF valence diagram (b) defined in ^ < a:: < 1 
region and the nonvalence diagram (c) defined in < x < ^ 
region. The large white and black blobs at the meson-quark 
vertices in (b) and (c) represent the ordinary LF wave function 
and the nonvalence wave function vertices, respectively. The 
small white blob at the quark-gauge boson vertex indicates 
the nonlocahty of the vertex. 

In Rcfs. [§,|lO|, the overlap representation of light-front 
wave function for the SQDs has been obtained from the 
formal definitions of light- front field operators. In this 
work, however, we shall derive the SQDs of the pion start- 
ing from the covariant Bethe-Salpeter (B-S) amplitude of 
the current J given by (see Fig. p|(a)) 



J" 



iNr 



d^k 5{x - k+/P+)H,o,H^^,S^' 
(27r)4 _ ^2 + ie][{k - A)2 - ^ 



[(P-fc)2 



(12) 



where Nc is the color factor, 5{x~k'^ / P^) represents the 
composite operator denoted by a small white blob at 
the quark-gauge boson vertex in Fig. ||, and i?cov {H'^^v) 
is the covariant initial (final) state meson-quark vertex 
function that satisfies the B-S equation. We refer to the 
black vertex in Fig. ||(c) as the nonvalence wave func- 
tion vertex as in the B-S formalism it is obtained by a 
continuation of the usual B-S amplitude, as we discuss 
below. 

Using the following identity 



+ m) + -j^l^{p~ 



Pon)> 



(13) 



we can express the trace term 5'' in Eq. (|T^) in terms of 
the on-mass shell propagating part of quark propagators 
and the instantaneous one as follows 



S^" = Tr 
= Tr 
+ Tr 

with 



75(1*1 + m)Y{A + "i)75(-^g + m) 
75(^10.1 + "^)7^(^2o„ + "^)75(-^gon + m) 



inst. 



(14) 




Tr 



inst. 



= 2(Pi -Pio 



2(P2 -P2on) 



M -t- -t- M 

P2onPqon ~ P2onPqon 



+5^^(P2on • Pqon + "1^) 



^'lon^^qon PlonPqon 



+5^^(Plon -Pqan+m^) 
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2(P„- -Pnon) 



M + I + M 
-PlonP2on + -PlonP2on 



-.9^^(Pion •P2on + rn^) 



+ V+Pfo„(Pr-Plo„)(P2 -P2on). (15) 

where pi = k,p2 = fc — A, and pq = P — k and the 
subscript (on) means on-mass shell quark propagator. 

The SQD is given by the solutions to the B-S equa- 
tion 

(Ml ~ Ml)x{x^,k,^) 

= [dy][(fl±]JC{x,,k,j_;yj,lj^)x{yj,lj±), (16) 



where K, is the B-S kernel which in principle includes all 
the higher Fock-state contributions, — M^/(l — ^), 
Ml = (m^ + k2_)/(l - x) - (m^ + k\)/{£_ - a;), and 
xixi,\ii±) is the B-S amplitude. Both the valence and 
nonvalence B-S amplitudes are solutions to Eq. (^6|). For 
the normal B-S amplitude, referred to as the valence wave 
function here, a; > ^, while for the nonvalence B-S ampli- 
tude X < ^. We use the notation for these two solutions 



X(2^2) X 
X(1^3) = X 



val 

nonval 



(17) 



This notation is motivated by the relationship to the 
Fock state picture, in which for the nonvalence vertex 
the parton number before and after the kernel is inter- 
preted as changing (from 1 to 3). However, as illustrated 
in Fig. |^(c), the nonvalence B-S amplitude is an analytic 
continuation of the valence B-S amplitude. In the LFQM 
the relationship between the B-S amplitudes in the two 
regions is given by [p^ 

(Ml - Ml)Xil-,3){x^,k^±) 

[dy][d^l±]/C(a;,, k,j_; yj, lji)X(2^2)(yj, IjJ-), (18) 



where again the kernel includes in principle all the higher 
Fock-state contributions because all the higher Fock com- 
ponents of the bound-state are ultimately related to the 
lowest Fock component with the use of kernel. This is 
illustrated in Fig. ||. 



Equations (|T^) and (^8|) are integral equations for 
which one needs nonperturbative QCD to obtain the ker- 
nel. We do not solve for the B-S amplitudes in this work, 
but a nice feature of Eq. (|lj) is a natural link between 
-^nonval ^vai ■^yjjjg}^ cnablcs an application of a light- 
front CQM even for the calculation of nonvalence contri- 
bution in Fig.j2(c). We describe this below. In 1 -I- 1 
QCD models |§19[, it is shown that expressions for the 
nonvalence vertex analogous to our form given in Eq. ( p^ ) 
are obtained. 

Corresponding to these two regions of the B-S ampli- 
tude, the Cauchy integration over k~ in Eq. (|l|) has 
two nonzero contribution to the residue calculations, 
one coming from the interval (I) A+ < /c+ < P"*" [see 
Fig. J|(b)] and the other from (II) < A:+ < A+ [see 
Fig-ic)]. 

(I) In the region of A+ < fc+ < P+, the residue is 
at the pole of fc" = P" - [k^^ + m'^ - ie]/{P - fc)+(i.e. 
Pq — PgoiJi which is placed in the upper half of complex- 
k~ plane. Thus, the Cauchy integration of J7+ in Eq. ( p^ ) 
over k~ gives 



^ dx 5{x-k+/P+) 



{P + P')+ 167r3 xil-x)x' 
X J t^^k_LX(2^2)(a;,ki)5+;X(2^2)(a;',k'L), 

(19) 



where 



X(2^2) 



Af2 - Af 2 



1-x 



^(2-2) - M2 - M^2 



and 



val ^-^ 



4P+ 



-(k^-k'^+TO^). 



(20) 



(21) 



The internal momenta of the (struck) quark for the final 
state are given by 



i-r 



k'± = kj 



1 



1-e 



(22) 




(l->3) "(2->2) 

FIG. 3. Nonvalence vertex(black blob) linked to an ordi- 
nary light-front wave function(white blob). 



While the light- front vertex function hi^p {h'^p) formally 
is given by the covariant i/cov (^cov); present 
work the radial wave function X(2— ►2) (X(2^2)) [conse- 
quently h^p {h'j^p)] is obtained from a light-front con- 
stituent quark model as we shall show later. It is also 
interesting to note that in this spectator pole diagram 
{Pq = Pqoa) with the plus component of the current, the 
instantaneous part in Eq. (15) does not contribute at all, 
i.e. all particles are on their mass shell. 

(II) In the region of < fc"*" < A"*", the residue is at the 
pole oik~ = [k^ + m2 — ie]/A;+ (i.e. p^ =Pion)' which is 
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placed in the lower half of complex-fc plane. Then the 
Cauchy integration of in Eq. ( |l2|) over fc^ reads 



« dx 5{x-k+/P+) 
{P + P')+ Jo 167r3 xil-x)x' 

dy 



2/(1 - y) 



d li/C(x,ki;y,U)x(2^2)(2/,l±), 



(23) 



where Eq. (|T^ has been used for the nonvalence wave 
function at the black blob in Fig. ||(c), and corre- 
sponds to the light-front energy denominator at the small 
white blob in Fig. ^(c). The explicit form of x^ is given 
by 



(1-0 



(24) 



where k" ± = kj_ + {x/£,)A±. W&,call x^ the hght-front 
vertex function of a gauge boson Id. In the calculation of 
the trace term 5*,+,, one can easily see from Eq. ( [l5| ) that 
there is one instantaneous contribution from the specta- 
tor (pg) line in addition to the on-mass shell contribution 
given by Eq. (^l|), which leads to 



9+ = 9^ 

i~' it-i /I) ^ „• 



AP 



(25) 




(b)+(c) 




FIG. 4. Effective treatment of the light-front nonvalence 
amplitude given by Fig. mc) . 



While one can in principle also consider the B-S amplitude 
for x^ : we note that such extension does not alter our results 
within our approximation in this work because both hadron 
and gauge boson should share the same kernel. 



Considering the quark-meson and quark-gauge boson 
vertices together, we also find that the four-body en- 
ergy denominator (-D4) appearing in Fig. |^ is absent in 
Eq. (^H). This absence of in our nonvalence cur- 
rent matrix amplitude given by Eq. ( p3| ) is due to the 
sum of two possible diagrams in the light-front time- 
ordering (see Figs, ^(b) and (c)). Summing over the 
two time-ordered diagrams Figs, lib) and (c), one can 
easily find the following identity, I/D4D2 + I/D4D2 = 
I/D2D2, which removes the complicate four-body en- 
ergy denominator term. We thus obtain the ampli- 
tude corresponding to the nonvalence contribution given 
by Eq. (23|) in terms of ordinary light-front wave func- 
tions of hadron (x(2^2)) and gauge boson (x^ ) as shown 
in Fig. §(d). This method, however, requires the rel- 
evant operator IC{x,'k±;y,l±) which is in general de- 
pendent on the involved momenta connecting the one- 
body to three-body sector as depicted in Fig. |3[ While 
the relevant operator IC is in general dependent on 
all internal momenta (x, kj_; y, lj_), the integral of K, 
over y and 1^ in Eq. (p3|), which we define as G^r = 
J[dy][d'^l±]K:{x,kj_;y,l±)x{2-*2){yA±), depends only on 
X and . In this work, we approximate as a constant 
which has been tested in our previous exclusive semilep- 
tonic decay processes |l^] and proved to be a good ap- 
proxiamtion at least for small momentum transfer region. 
As we shall show in the next section, the validity of this 
approximation can be checked by examining the frame- 
independence of our numerical results. 



IV. MODEL CALCULATION OF SQD 

In the previous sections we have derived the skewed 
quark distribution function x, t) starting from a co- 

variant model. In this section we replace the light-front 
vertex function h^p (or equivalently X(2^2xL^h^i-)) by 
the standard light-front vertex function |ri|,|l2|,Q , which 
is symmetric in the variables of the constituent qq pair 
and has been successful in predicting many static prop- 
erties of ground state mesons. Different choices of the 
vertex function are of course possible. 

Comparing Xi;2^2) with our light-front wave function 
given by Ref. [p2[ , we identify 




X{2^2){x,k^) 



where the Jacobian of the variable tranformation k = 
{kz,k±) — > (a;,kj_) is obtained as dkz/dx — Afo/[4a;(l — 
x)] and the radial wave function is given by 



^(k^) = 



1 



■exp(-kV2/32), 



^3/2^3 

which is normalized as / (i^fc|(/)(k^)p — 1. 



(27) 
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Substituting Eqs. (gej) and (|27|) into Eqs. (|19|) and (|23|), 
we obtain the valence and nonvalence contributions to the 
SQDs of the pion in LFQM 



1-1 



(a;,k^)- 



dx' 



dx 



0(x',k\) 
. (28) 



and 



^r(C,^,i) - / rf^k^^^</.(x,k^)x^(x,kV) 

kj^ • k'i + m2 + a;(l - a;)a;'(M2 - M^) 



a;'(l-a;')V'k^ 



(29) 

where we treat the last term in Eq. (|2^) as a constant G^r 
which will be fixed by the sum rule expressed in terms of 
Tl''^ and as 



j\xTr\^,x,t)+ J^xT!," 



i^,x,t), (30) 



for given —t. We note that Eq. ( |30| ) is used as a constraint 
on the frame-independence of our model. 

In our numerical calculations, we use the model param- 
eters (m, /?) = (0.22, 0.3659) [GeV] obtained in Rcf. 
for the linear confining potential model. Before we cal- 
culate the SQDs and the form factor of the pion, we 
first consider the valence contribution to the pion EM 
form factor and see how much the nonvalence contribu- 
tion is needed to obtain the frame-independent result of 
our model. To this end, we show in Fig. ^ the valence 
contribution (see Fig. g(b)) to the pion EM form fac- 
tor with different values of ^, where the thick solid, cross 
(x), dotted, dot-dashed, and thin solid lines represent the 
purely transverse (^ = and ^ 0), ^ = 0.02, 0.3, 0.6, 
and the purely longitudinal (i.e. ^ 7^ and A± — 0) 
result s, resp ectively, and compare with the experimental 
data [g2 23|. Note that the purely transverse frame result 
(thick solid line) is the exact solution within our model 
calculation. As one can see from Fig. ^, the nonvalence 
contribution to the pion form factor, i.e. the difference 
between ^ = and ^ ^ results, increases as ^ does. 
Our special interesting region in this work is the small 
—t region where the nonvalence contribution is especially 
large and our effective method for the calculation of the 
nonvalence contributions works pretty well. 

Including the nonvalence contribution given by 
Eq. ( p9[ ) to the SQDs of the pion, we can determine 
our constant G^r in a frame-independent way by the 
sum rule dx J^™'(C = 0,x,t) = dx J^r'(C 



0,x,t) + J^dx ^r(e ^ 0,x,t) = F^{t). In Fig. I, we 
show the ^-dependence of Gtt for different — t-values, i.e. 
—t = (diamond), 0.2 (black circle), 0.5 (white circle), 
and 1.0 (black square) [GeV^], respectively. As one can 
see in Fig. ^, G-^ shows approximately constant behavior 
for ^ > 0.1 at given small —t. It is not surprising to see 
that Gtt becomes very large as ^ ^ 0, however, this does 
not cause a significant error in our calculation because 
the nonvalence contribution in the very small ^ region is 
highly suppressed. Note from Eq. ( p3[ ) that J^™ has the 
form of = G^ X thus Gjr must be very large as 

^ — + to give the small contribution of , since the in- 
tegral vanishes. Therefore the results arc consistent with 
an almost constant value for G^r at least for small —t as 
we show below. On the other hand, we note that there is 
an obvious i-dependence for G^r , which might be the limit 
of our constant approximation. In principle, we can ob- 
tain the SQDs in a frame-independent way by using the 
true values of Gjr as shown in Fig. || for given (^,t). In 
the following, we compare the SQDs and the form factor 
obtained from true values of Gjr (i.e. frame-independent 
result of our model) with those obtained from a single 
average value of G^r — Gave. = 0.32 for all (^, t) to check 
the reliability of our constant G,r approximation. 

In Figs. and ||, we show the SQDs J>r(^,x,t) of 
the pion for fixed momentum transfer —t = 0.2 GeV^ 
(0 < C < 0.92) and -t = 1.0 GeV^ (0 < ^ < 0.98) 
but with different skewedness parameters f , respectively. 
The solid and cross (x) lines in the nonvalence contribu- 
tions are the exact solutions obtained from true values 
of Gtt and the effective ones obtained from our average 
value of Gave. = 0.32, respectively. The dotted lines rep- 
resent the instantaneous contributions to JP^^(^, x, t) ob- 
tained from true values of G^r. The SQDs at ^ = as 
shown in Figs. 0(a) and ^(a) correspond to the ordinary 
quark distributions with vanishing nonvalence contribu- 
tions. The frame-independence of our model calculation 
is ensured by the area under the solid lines (valence -I- 
nonvalence) being equal to the pion form factor at given 
—t. As one can see from Figs. |7|(b-c) and ||(b-c), while 
the nonvalence contributions are small for small f ~ 0.3, 
they are large for large skewdness parameter — 0.9. 
In our model calculations, the nonvalence contributions 
obtained from true values of G^r (solid lines in each fig- 
ure) at -t=0.2 (1.0) GeV^ for ^=0.3 (0.3) and 0.9 (0.9) 
are approximately 11 (4) % and 90 (85) %, respectively. 
Comparing with the exact solutions, the numerical re- 
sults with a single average Gave. = 0.32 (cross lines in 
each figure) are shown to reproduce the exact ones up 
to 97 % for ^ = 0.3 and 90 % for ^ = 0.9, respectively. 
It is also interesting to note that the instantaneous con- 
tribuitons (dotted lines in each figure) become more pro- 
nounced as ^ ^ Cmax for each —t, which is a very dif- 
ferent feature from a scalar theory model where there 
is no such instantaneous contribution. While the instan- 
taneous part of the nonvalence contribution vanishes as 
X — > = lim^^oiS,^^) as shown in Figs, ^and^, the net 
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result of x, t) Jncluding the on-mass shell propa- 

gating part does notEl vanish as a; = ^ and consequently 
causes a discontinuity to the zero value of 
However, such discontinuity at x = ^ is just an artifact 
due to the difference in the x ^ behavior between 
the gauge boson vertex (^^(a;,]^")) in Eq. and the 
hadronic vertex {x2'^2i^' ^ Eq. (19) of our approx- 

imate model calculation. We have indeed confirmed that 
the discontinuity at x = ^ does not occur in the limit of 
a point hadcon vertex as already noticed in the QED cal- 
culation l^u. Thus, for a full analysis of DVCS satisfying 
the factorization theorems , it would be necessary to 
solve the bound-state B-S equation similar to Eq. (|l6|) for 
the gauge boson (x^) as well as for the hadron (x2^2)- 
In the chiral quark-soliton model analysis of the nucleon 
SQDs 1^^, the discontinuity at a; = ^ was imputed to 
an artifact of neglecting the momentum dependence of 
the constituent quark mass. Both the B-S amplitude 
for and the dynamical quark mass would be anyway 
necessary for the model improvement. Nevertheless, our 
effective method seems useful for the present study of 
the relation between SQDs and the form factor in the 
nonperturbative regions. 

In Fig. ^, we show our effective solutions of the pion 
form factor for ^ = 0.3 (thin solid line) and 0.9 (long- 
dashed line) cases obtained from our average value of 
Gtt = Gave. = 0.32 and compare with the exact solution 
with ^ = (thick solid line) as well as the experimental 
data [^,^. The dotted and dot-dashed lines represent 
the valence and the instantaneous contributions to the 
form factor for the case of ^ = 0.9, respectively. In fact, 
there are — tmin values for nonzero ^ due to A^^ > 
(see Eq. (H)). We thus use the analytic continuation by 
changing A_l to iA± in Eqs. ( p^ ) and ( |29| ) to obtain the 
result for < —t < — imin where there is no singular- 
ity. A continuous behavior of the form factor near — imin 
confirms the analyticity of our model calculation. Our 
effective solution (thin solid line) with ^ — 0.9 shows al- 
most maximum deviation (< 10 %) from the exact one 
(thick solid line) and the deviation becomes smaller as 
^ reduces. Our effective method of evaluating the non- 
valence diagram with a constant operator Gtt shows a 
definite improvement to restore the frame-independence 
of our model and seems to be a quite reliable approxima- 
tion. 

In Fig. [l^, we show the nth moments(n = 1,2,3) of 
^7r(Cj ^1 given by Eq. ( pT| ) a,t t — using the true value 
of Gtt shown in Fig. ^. Although at t = the skewedness 
parameter vanishes (see Eq. (|^)), we use the model forms 
for Jf"™' and J^;'" given by Eqs. (|2|,|2|) to define the ex- 
trapolation to ^ ^ 0. The thick solid, dotted, and dot- 
dashed lines represent the total(=valence -|- nonvalence), 



valence, and nonvalence contributions, respectively. For 
comparison, we also show the nonvalence contribution 
obtained from the average value of G-^ — 0.32(cross 
lines). As one can see, the first(n — 1) moment (top 
thick solid line) is ^-independent because the sum rule 
for n = 1 yields the physical pion form factor, Fi (^, t) — 
FTr{t). Also, the higher moments F2(^, i) (middle thick 
solid line) and F3(^, i) (bottom thick solid line) satisfy 
the polynomiality conditions(See Eq. (pT|)) discussed in 
the previous section. In Fig. |l^ we also plot the phe- 
nomenological form 



^ Fni^ ^ 1) - (1 - O'lFniC 



l)-i^n(e-O)], (31) 



shown by the diamonds. The numerical values for i^^ = 
Fn{£, — 0) and F^"" — F„{£, ~ 1) are summarized in Ta- 
ble I. Our value of F^ is in good agreement with the value 
obtained by the QCD sum rules with non-local conden- 
sates m. 

We also compute the nth moment of the tt wave func- 
tion, 0,(2/) = ^^(^ = 0,y,i = 0), defined by pi 



dyy^'My) 



(32) 



^In fact, thisbehavior of J-^'" {^,x,t) at x = ^ has been 
anticipated in |9| without a proof. 
^We are grateful to M. Diehl for discussion of this point. 



where y — 2x — 1. Our numerical results are summarized 
in Table II and compared with several other theoretical 
results. It is interesting to note that our value of (1/2) 
is very close to the asymptotic value (0.2) of the sec- 
ond moment obtained from the well-known asymptotic 
quark distribution amplitude 4>-K{y) = f (l~y^)i which is 
quite different from the early work obtained with QCD 
sum rules shown in the second row of Table II. A 
later calculation using QCD sum rules p7| , shown in the 
third row, found moments closer to the asymptotic val- 
ues. However, we emphasize that the Bethe-Salpeter am- 
plitude used in the present work is a model rather than a 
solution to the B-S equation with a kernel derived from 
nonperturbative QCD. 

Finally, in Fig, 11, we obtain the isosinglet SQDs of the 
pion by subtracting the valence part ^7r(^ = 0, a:;,i = 0) 
from the nonvalence part .^^(^ = l,x,i = 0). Our result 
(dotted line) in Fig. |ll| is qualitatively very similar to 
the total isoscalar skewed quark distribution in the pion 
satisfying the soft pion theorem (Fig. 5 of Ref. |^), that 
is obtained by the low-energy effective field theory based 
on the instanton model of the QCD vacuum. 



V. CONCLUSION 

In this work, we investigated the SQDs of the pion for 
small momentum transfer (— t < 1 GeV^) region in the 
light-front quark model. Since the light-front nonvalence 
contributions to the SQDs of the pion are large especially 
at small momentum transfer region as shown in Fig. 
it is very crucial to take them into account to guarantee 
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the frame-independence of the model. Applying our ef- 
fective treatment [|l6| , i.e. the nonvalence B-S amplitude 
given by Eq. of the nonvalence contribution to the 
SQDs JSr {^r^ + TJ^^') of the pion, we express T^" (see 
Eqs. (^3|) and (p9|)) in terms of ordinary light-front wave 
functions of a gauge boson and a hadron and calculate 
this nonvalence contribution numerically. 

The main approximation in our effective calculation 
is the treatment of relevant operator A^(a;, k^; y, Ij^) in 
Eq. ( |2^ ) connecting the one-body to three-body sec- 
tor (see Fig. ^) by taking a constant via G^r = 
/[dy][(i2l_L]/C(x,k_L;y,l_L)x(2^2)(y,l_L), which in general 
depends on x and k^. The reliability of this constant 
approximation was checked by examining the frame- 
independence of our numerical results using the sum rule 
given by Eq. (|30|), i.e. the exact results of TTr{^,x,t) 
and F7r(t) obtained from the true values of G^r given by 
Fig. ^ were compared with those obtained from our sin- 
gle averge value of G^ = 0.32 for all (^,t). The numeri- 
cal results of our constant G^ prescription have shown 
definite improvement (better than 90 % accuracy for 
^ < 0.9) to restore the frame-independence of our model 
(see Figs. 0- ^ and seemed to be a quite reliable ap- 
proximation. Our model also satisfies the polynomiality 
conditions for n = 2 and 3 moments of the SQDs (See 
Fig. and our value of is in good agreement with 
the QCD sum-rule result with non-local condensates [ p5| . 
Moreover, the ordinary quark distribution amplitude of 
the pion ((/'^(y — 22; — 1)) can be obtained from the 
two-pion distribution amplitude in ^ = and t — limit 
where one of the produced pions becomes soft. Our result 
of (pTT is very close to the asymptotic quark distribution 
amplitude consistent with the CLEO measurements . 
The isosinglet skewed quark distribution amplitude in 
the pion (See Fig. ^l]) is also consistent with the result 
satisfying the low-energy soft pion theorem ||^ . 

For the model improvement, however, it would be nec- 
essary to consider not only the dynamical quark mass 
discussed in our previous work but also the B-S am- 
plitude for the quark-gauge boson vertex to remove the 
apparent discontinuity at x = ^ as shown in Figs.^ and||. 
Especially, in exploring the new physics associated with 
the SQDs for the treatment of the pion form factor, we 
note from Eq. (^ that the definition of J-V involves the 
matrix element of the light-front operator present in the 
three-point approach to quark distributions of hadrons in 
scaling regions. As was shown in Ref. ||3^ the treatment 
of such operators can lead one to consider nonlocal quark 
condensates, which can introduce nonperturbative QCD 
structure in the quark-gauge boson vertex. This will be 
a subject of future research by the authors. 
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FIG. 5. The valence contribution to the pion EM form fac- 
tor with different skewedness parameters ^ compared with the 
experimental data p2lJ2l3|. 
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FIG. 6. The ^-dependence of G-n for different momentum 
transfers —t = (diamond), 0.2 (black circle), 0.5(white cir- 
cle), and 1 (black square) [GeV^], respectively. 



FIG. 7. Skewed quark distributions of the pion at —t = 0.2 
GeV^ with ^ = in (a) ,0.3 in (b), and 0.9 in (c), respectively. 
The solid [cross (x)] line in nonvalence contribution represents 
the full result of using true [average] G,r value and the dot- 
ted line represents the instantaneous part of the nonvalence 
contribution. 
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FIG. 9. The effective solution of pion form factor using a 
common average G-^ = C?ave. = 0.32 value for ^ = 0.3 (tliin 
solid line) 0.9 (long-dashed line) compared with the exact 
solution (thick solid line) as well as experimental data [ p^ , p3| . 
The dotted and dot-dashed lines represent the valence and 
instantaneous contributions to the form factor for ^ = 0.9 
case, respectively. 
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FIG. 8. Skewed quark distributions of the pion at — t = 1 
GeV^ with different ^ = in (a), 0.3 in (b), and 0.9 in (c), 
respectively. The same line code is used as in Fig. ^ 
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FIG. 10. The nth moments F„{^,-t) of SQDs of the 
pion at t = using the true value of G. The total re- 
sults(thick solid lines) for n = 2 and 3 are well fitted by 
the simple polynomial(diamond) in 5. For comparison, we in- 
clude the nonvalence contributions obtained from the average 
Gtt — 0.32(cross lines). 
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FIG. 11. The SQDs of the pion at t = using 
the true value of G(=0.384). The thick solid, solid, 
and dotted lines represent J-tt{£, ~ 0,x,t = 0) (valence), 
^7r(^ = l,x,t — 0) (nonvalence), and their difference 
■^vr(^ = l,x,t = 0) ~ ^wi^ = 0,x,t = 0), respectively. 



TABLE I. The nth moments Fn{^, x,t = 0) of the SQDs for the pion where F,X = F„ {^ = 0) and 
F^^ = Fn{(, — 1) represent pure valence and pure nonvalence contributions, respectively. 



Model 


F^Fn 


F-nF-n F^Fn F^iFn F^iFn 


F:^[Fr] 


( 

I 


Durs 
^ 


0.503[0.623] 


0.312[0.433] 0.215(0.322] 0.16[0.25] 0.123(0.201] 
0.25[-] 
0.29H 


0.098(0.165] 


TABLE II. The nth moments (y„) of the ordinary quark distribution amplitude 
pion. 


^,r(y) for the 


Model 


(yi) 


(^2) (y:i) (yi) (ys) 




( 

1 
1 
1 


3u 
26 
27 

28 


rs 

C 


0.0 
0.0 
0.0 
0.0 


0.239 0.0 0.109 0.0 
0.43 0.0 0.24 0.0 
0.25 0.0 0.12 0.0 
0.25 0.0 0.11 0.0 


0.062 
0.15 
0.07 
0.06 



For the gaussian parameter /3 — 0.36 GeV in Ref. p,8\. 
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